Divergent Time Scale in Axelrod Model Dynamics 
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We study the evolution of the Axelrod model for cultural diversity. We consider a simple version 
of the model in which each individual is characterized by two features, each of which can assume q 
possibilities. Within a mean-field description, we find a transition at a critical value q c between an 
active state of diversity and a frozen state. For q just below q c , the density of active links between 
interaction partners is non-monotonic in time and the asymptotic approach to the steady state is 
controlled by a time scale that diverges as (q — (fc) -1 ^ 2 ■ 
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A basic feature of many societies is the tendency to 
form distinct cultural domains even though individuals 
may rationally try to reach agreement with acquain- 
tances. The Axelrod model provides a simple yet rich 
description for this dichotomy by incorporating societal 
diversity and the tendency toward consensus by local in- 
teractions . In this model, each individual carries a set 
of F characteristic features that can assume q distinct 
values; for example, one's preferences for sports, for mu- 
sic, for food, etc. In an elemental update step, a pair of 
interacting agents i and j is selected. If the agents do 
not agree on any feature, then there is no interaction. 
However, if the agents agree on at least one feature, then 
another random feature is selected and one of the agents 
changes its preference for this feature to agree with that 
of its interaction partner. A similar philosophy of al- 
lowing interactions only between sufficiently compatible 
individuals underlies related systems, such as bounded 
confidence and constrained voter- like models Q . 

Depending on the two parameters F and q, a phase 
transition occurs between cultural homogeneity where 
all agents are in the same state, and diversity 0, 01 S 
6]. The latter state could either be frozen, where no 
pair of interacting agents shares any common feature, 
or it could be continuously evolving if pairs with shared 
features persist. The rich dynamics of the model does 
not fall within the classical paradigms of coarsening in 
an interacting spin system or diffusive approach to 
consensus in the voter model ||. In this Letter, we solve 
mean-field master equations for Axelrod model dynamics 
and show that the approach to the steady state is non- 
monotonic and extremely slow, with a characteristic time 
scale that diverges as q — > q c (Figs.QJ&|2J. 

The emergence of an anomalously long time scale is un- 
expected because the underlying master equations have 
rates that are of the order of one. Another important 
example of wide time-scale separation occurs in HIV [9j . 
After an individual contracts the disease, there is a nor- 
mal immune response over a time scale of months, fol- 



lowed by a latency period that can last beyond 10 years, 
during which an individual's T-cell level slowly decreases 
with time. Finally, after the T-cell level falls below a 
threshold value, there is a final fatal phase that lasts 2-3 
years. Our results for the Axelrod model may provide 
a hint toward understanding how widely separated time 
scales arise in these types of complex dynamical systems. 
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FIG. 1: Master-equation time dependence of bond densities 
Po, Pi, and P2 for q — q c — 4 _1 . Each agent has 4 neighbors. 
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FIG. 2: Master-equation result for Pi(t) for q = q c — ^T k , with 
k = —1,1,3,5, and 7 (progressively lower minima). Each 
agent has 4 neighbors. The dashed line has slope —2 (see 
text). Inset: Same data on a linear scale with T = t(q — q c ) 1 ^' 2 . 
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Following Refs. , we describe the Axelrod model in 
a minimalist way by the density P m of bonds of type to. 
These are bonds between interaction partners in which 
there are to common features. This description is conve- 
nient for monitoring the activity level in the system and 
has the advantage of being analytically tractable. We 
consider a mean-field system in which each agent can in- 
teract with a fixed number of randomly-selected agents. 
Agents can thus be viewed as existing on the nodes of 
a degree-regular random graph. Such a topology is an 
appropriate setting for cultural interaction, where both 
geographically nearby and distant individuals may in- 
teract with equal facility. We verified that simulations 
of the Axelrod model on degree-regular random graphs 
qualitatively agree with our analytical predictions, and 
this agreement becomes progressively more accurate as 
the number of neighbors increases (Fig. UJ. Thus the 
master equation approach describes the Axelrod model 
when random connections between agents exist. 
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FIG. 3: Active bond density from the master equations 
(curves) and from simulations of 10 2 realizations (A) on a 
degree-regular random graph with 10 4 nodes for various co- 
ordination numbers, and q = 8 states per feature. Inset: One 
realization with coordination number 4. 

If interaction partners share no common features (m = 
0) or if all features are common (to = F), then no inter- 
action occurs across the intervening bond. Otherwise, 
two agents that are connected by an active bond of type 
m (with < to < F) interact with probability m/F, af- 
ter which the bond necessarily becomes type m + 1. In 
addition, when an agent changes a preference, the index 
of all bonds attached to this agent may either increase 
or decrease (Fig. UJ. The competition between these di- 
rect and indirect interaction channels underlies the rich 
dynamics of the Axelrod model. 

Because we obtain qualitatively similar behavior for 
the density of active links, P a = J2k=i ^k, for all F > 2, 
we focus on the simplest non-trivial case of F = 2. For 
this example, there are three types of bonds: bonds 
of type (no shared features) and type 2 (all features 
shared) are inert, while bonds of type 1 are active. As 
q —* q c from below, P± is non-monotonic, with an in- 
creasingly deep minimum (Fig. l2*jl . while for q > q c , Pi 
decays to zero exponentially with time. There is a dis- 
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FIG. 4: Illustration of the state-changing bond updates when 
agent j changes state from ai&2 — > a\b\. The values at the 
right give the relative rates of each type of event. 



continuous transition at q c from a stationary phase where 
the steady-state density of active links P* is greater than 
zero to a frozen phase where Pf = 0. 

When fluctuations are neglected, the evolution of the 
bond densities P m when a single agent changes its state 
is described by the master equations: 



dPo 
dt 
dPi 
~dt 
dP 2 
~~dt 



V 



?7+l 
Pi 



Pi 



1 



-AP + 



n 



r)+l r)+l 
Pi V 



Pi 



rj+1 rj+1 



Pi 



AP - ^Pi + P 2 



^P-P 2 
2 



(1) 
(2) 
(3) 



where 77 + 1 is the network coordination number. The 
first term on the right-hand sides of Eqs. J2J and @ 
account for the direct interaction between agents i and j 
that changes a bond of type 1 to type 2. For example, 
in the equation for =Q, a type-1 bond and the shared 
feature across this bond is chosen with probability Pi/2 
in an update event. This update decrements the number 
of type-1 bonds by one in a time dt = -k, where N is the 
total number of sites in the system. Assembling these 



factors gives the term 



Pi 

ri+l 



Eq. ©. 



The remaining terms in the master equations represent 
indirect interactions. For example, if agent j changes 
from (ai, 62) to (di, 61) then the bond to agent k in state 
(ai,6i) changes from type 1 to type 2 (Fig. 0J. The 
probability for this event is proportional to PiA/2: Pi 
accounts for the probability that the indirect bond is of 
type 1, the factor 1/2 accounts for the fact that only the 
first feature of agents j and k can be shared, while A is 
the conditional probability that i and k share one feature 
that is simultaneously not shared with j. If the distribu- 
tion of preferences is uniform, then A = (q— l) -1 . As the 
system evolves A generally depends on the densities P m . 
Here we make an assumption of a mean-field spirit that 
A stays constant during the dynamics 5] ; this makes the 
master equations tractable. Our simulations for random 
graphs with large coordination number match the master 
equation predictions and give A nearly constant and close 
to (q — 1) _1 (Fig.[SJ, thus justifying the assumption. 

Let us first determine the stationary solutions of the 
master equations. A trivial steady state is Pf = 0, cor- 
responding to a static society. A more interesting sta- 
tionary solution is Pf > 0, corresponding to continuous 
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evolution; as we shall see, this dynamic state arises when 
q < q c . Setting = in the master equations and 
solving, we obtain: 
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Since A = X(q) is the only parameter in the master 
equations, the two stationary solutions suggest that there 
is a transition at a critical value q c such that both solu- 
tions apply, but on different sides of the transition. To 
locate the transition, it proves useful to relate Pi and P 2 
directly. Thus we divide Eq. (J3J) by Eq. and eliminate 
Pq via Pq = 1 — Pi — P2 and obtain, after some algebra: 



dPi __ -I + 77A- ±77(1 
dP 2 ~ " 



3A)Pi + 77(1 - A)P 2 



1 



jr?APi - t/P, 
The solution to Eq. JSJ) has the form 



Pi = a + j3P 2 - ^7 + SP 2 , 



(5) 



(G) 



where we determine the coefficients a, (3, 7 and S by 
matching terms of the same order in Eq. (JSJ) and in 
from Eq. JBJ). The procedure gives the solution except 
for one constant that is specified by the initial condi- 
tions. For the initial condition where features for each 
agent are chosen uniformly from the integers [0, q — 1], 
the distribution of initial bond densities is binomial, 
p m(t = 0) = ^!jlW( 1 /9) m (l - l/q?- m - Matching 
this initial condition to the solution of Eq. © gives: 
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As a function of P 2 , Pi has a minimum P™ m {q) that 
monotonically decreases as q increases and becomes neg- 
ative for q larger than a critical value q c (Fig. [5J| . The 
phase transition between the active and the frozen state 
corresponds to the value of q where Pi first reaches zero. 
To find q c , we calculate p™ m as a function of X(q) from 
Eq. (0) and then find the value of q at which p™ m be- 
comes zero. This leads to 



pmin 



47 ? A-(1 + A) 2 _ S(\,t}) 



2r/(l + A) 2 277(1 + A) 2 
from which the critical point is given by 



q c = 277 + 2^r){ii - 1) , 

while P™' 11 oc S oc (q c - q) for q < q c . 

We now determine the steady-state bond densities in 
the frozen state. From Eq. {7|), we compute the station- 
ary value Pf at the point where Pi first reaches zero. 
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FIG. 5: Pi vs P 2 from Eq. for 77 = 3 and q = q c - 2, q c - 1, 
9c, q c + 1 and q c + 2 (top to bottom). 



The smallest root of this equation then gives 
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1 + A + 277A - y/(l + A) 2 - 477A 



27/(1 + A) 



while P S = 1 - P|. 



The most interesting behavior is the time dependence 
of the density of active bonds, Pi(t). We solve for Pi(t) 
by first inverting Eq. (JJJ to express P 2 in terms of Pi 
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and then writing Pq = 1 — Pi — P 2 (Pi) also in terms of 
Pi , and finally substituting these results into the master 
equation J2J for Pi . After some algebra, we obtain 



dPi 



= SPi - (1 - A)v/277(l + A) 2 Pi-5 Pi 
- 27/(1 + A) 2 P 2 , 



(8) 



where we use the rescaled time variable t = 2 (ri+i){i+x) • 
This master equation can be simplified by substituting 
the quantity A = 277(1 + A) 2 Pi — S, which measures the 
deviation of Pi from its minimum value, in Eq. JSJ). We 
obtain 



dA 
dr 



-\/A(5 + A)(l-A + v / A) 



(9) 



Performing this integral by partial fraction expansion 
gives 




f S + A 

Ua(i-a) 2 

(V=ff-1-A)( 



- 2 In 1 ± 




(10) 



For q > q c , only the upper sign is needed. For q < q c , 
the upper sign applies for t < i mm and the lower sign 
applies for t > /j min ; here i min is the time at which P\(t) 
reaches its minimum value. Substituting back t and Pi 
in Eq. (|10|l gives the formal exact solution of Eq. JHJl . 
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For q < q c , we determine P\{t) near its minimum by 
taking the A —> limit of Eq. @ . This gives 



dA 
~dt 



iSVA, 



(11) 



with a 



(i-A) 



> 0. For S > 0, the solution 



2(r,+l)(l + A) 

to the lowest-order approximation shows that Pi has a 
quadratic form around its minimum: 



Pi(t) - 



- P™ in oc A ; 



' 877(1 



^^(t-W) 2 . (12) 



When q — > <7 C , the factor 5 may be neglected as long as 
A > S, and this leads to A decaying as t~ 2 before the 
minimum in Pi is reached (dashed line in Fig.__|). 

The peculiar behavior of Pi as a function of time for 
q below but close to the critical value q c is shown in 
Fig. [__ The density of active bonds quickly decreases 
with time and this decrease extends over a wide range 
when q is close to q c . Thus on a linear scale, Pi remains 
close to zero for a substantial time. After a minimum 
value at i m i_ is reached, Pi then increases and ultimately 
reaches a non-zero asymptotic value for q < q c . The 
quasi-stationary regime where Pi remains small is defined 
by: (i) a time scale of the order of one that characterizes 
the initial decay of Pi{t), and (ii) a much longer time 



scale t a 



, where Pi rapidly increases and then saturates 



at its steady-state value. 

We can give a partial explanation for the time depen- 
dence of Pi. For q > q c , there are initially small enclaves 
of interacting agents in a frozen discordant background. 
Once these enclaves reach local consensus, they are in- 
compatible with the background and the system freezes. 
For q < q c (less diversity), sufficient active interfaces are 
present to slowly and partially coarsen the system into 
tortuous domains whose occupants are either compati- 
ble (that is, interacting) or identical. Within a domain 
of interacting agents, the active interface can ultimately 
migrate to the domain boundary and facilitate merging 
with other domains; this corresponds to the sharp drop in 



Po seen in Fig. [T| 01 ■ While this picture is presented in 
the context of a lattice system, it remarkably still seems 
to apply for degree-regular random graphs and in a mean- 
field description. 

Both i m ; n as well as the end time of the quasi- 
stationary period t a symp increase continuously and di- 
verge as q approaches q c from below. To find these diver- 
gences, we expand t mm and t a symp in powers of S. From 
Eq. i|10|) , the first two terms in the expansion of t n 
5^0, are 



as 



t min = t(Pr n )^AlnS- 
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as the time at 



(Qc-q) 



where A, B are constants. As a result, t„ 
as q — > q c . Similarly, we estimate i aS y m 
which Pi reaches one-half of its steady-state value. Using 
Eqs. (gj and t_Tj|), we find 
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as q 



as S , 
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For q > q c , the system evolves to a frozen state with 
Pi -> 0. To lowest order Eq. © becomes ^ = 



dt 

Here T > since S < for 



with T - 2(t?+1)(1 + A) 

q > q c . Consequently Pi decays exponentially in time as 
t — > 00. As q approaches q c , S asymptotically vanishes 



as (q c — q) and the leading behavior is T ~ (q — q c 



-1/2 _ 

Thus again there is an extremely slow approach to the 
asymptotic state as q approaches q c . 

In summary, the density of active links is non- 
monotonic in time and is governed by an anomalously 
long time scale in the 2-feature and q preferences per fea- 
ture Axelrod model. For q < q c , an active steady-state 
state is reached in a time that diverges as {q c — q)^ 1 ^ 2 
when q —> q c from below. For q > q c , the final state is 
static and the time scale to reach this state also diverges 
as (q c — q)^ 1 / 2 as q — > q c from above. 
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